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Explicit expressions for the non-singlet and singlet spin-dependent structure function g 1 in the 
small- x region are obtained. They include the total resummation of the double- and single- loga- 
rithms of x and account for the running QCD coupling ct s effects. Both the non-singlet and singlet 
structure functions are Regge behaved asymptotically, with the intercepts predicted in agreement 
with experiments. A detailed comparison with the DGLAP evolution equations for different values 
of x and Q 2 is performed. Finally, the role played by singular terms in DGLAP fits for the initial 
quark densities is discussed and explicitly shown to mimic the resummation of leading logarithms 
at small- a;. 
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I. INTRODUCTION 
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In the standard theoretical approach for investigating the DIS structure function gi(x,Q 2 ), namely DGLAP 0, 
gDGLAP j g a convolution of the coefficient functions Cdglap and the evolved parton distributions. The latter are 
also expressed as a convolution of the splitting functions Pdglap and initial parton densities, which are fitted from 
experimental data at large x, x ~ 1 and Q 2 ~ 1 GeV 2 . However there is an obvious asymmetry in treating the 
• Q 2 - and x- logarithmic contributions in DGLAP. Indeed, the leading Q 2 - contributions, ln(Q 2 ), are accounted to all 
' orders in a s whereas Cdglap(x) and Pdglap(x) are known in first two orders of the perturbative QCD. On the 
Oh! other hand, in the small- a: region the situation looks opposite: logarithms of x, namely double logarithmis(DL), i.e. 
O ,' the terms (a s \n 2 (l/x)) k , and single logarithms (SL), the terms (a s ln(l/x)) fc ,with k — 1,2,.., are becoming quite 
sizable and should be accounted to all ordersi in a s . The total resummation of DL terms was first done[2| in the 
fixed a s approximation, and led to a new expression g® L , for g±, that in the small- x asymptotics was of the Regge 
(power-like) form and much greater than the well-known small- x asymptotics of g® GLAP . 

Strictly speaking, the results of Refs. could not be compared in a straightforward way with DGLAP because 
instead of the running a s , with the parametrization 



DGLAP (rfl\ f1\ 

a s = a s {Q ), (1) 

Refs. |2( had used a s fixed at an unknown scale. A closer investigation of this matter|5j led us to conclude that 
the DGLAP- parametrization of Eq. (yi can be a good approximation at x not far from 1 only. Instead, a new 
parametrization was suggested, where the argument of a. in each of the ladder rungs of the involved Feynman 
diagram is the virtuality of the horizontal gluon (see Ref. p| for detail) . Indeed this parametrization works well both 
for small and large x. It converges to the DGLAP- parametrization at large x but differs from it at small x, and it 
allowed us to obtain in Refs. [?J the expressions for g\ accounting for all-order resuummations of DL and SL terms, 
including the running a s effects 1 . This led us to predict the numerical values of the intercepts of the singlet and 
non-singlet g\. These results were then confirmed @ by several independent groups who have analyzed the HERMES 
data and extrapolated them at small x. 

On the other hand, it is well known that, despite missing the total resummation of In a;, DGLAP works quite 
successfully at x <C 1. This might suggest that the total resummation of DL contributions performed in Refs. 
should not be relevant at available values of x and might be of some importance at extremely small x reachable in the 
future. In Ref. |9( we made a detailed numerical analysis and explained why DGLAP fits can be successful at small 



1 The parametrization of see Ref. |^ was used later in Refs. [f| for studying the small-a: contribution to the Bjorken sum rule. 
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x. Indeed in order to describe the available experimental data, singular expressions (see for example Refs. [ToL ITU') 
are introduced for the initial parton densities. These singular factors (i.e. the factors which — > oo when x — ■> ) 
introduced in the fits mimic the total resummaton of Refs. Q. Then using the results of Ref Q for incorporating the 
total resummation of In x it allows to simplify the rather sophisticated structure of the standard DGLAP fits down 
to a normalization constant at small x. 

II. DIFFERENCE BETWEEN DGLAP AND OUR APPROACH 

In DGLAP, gi is expressed through convolutions of the coefficient functions and evolved parton distributions. As 
convolutions look simpler in terms of integral transforms, it is convenient to represent g\ in the form of the Mellin 
integral. For example, the non-singlet component of g\ can be represented as follows: 



<h S DGLAp(x,Q 2 ) = (eJ/2) / 3 ^{l/xyC DGLAP {uj)8q{uj)exv 



dk 2 



(2) 



with Cdglap{u) being the non-singlet coefficient functions, Jdglap{u, a s ) the non-singlet anomalous dimensions 
and 5q{uj) the initial non-singlet quark densities in the Mellin (momentum) space. The expression for the singlet g\ 
is similar, though more involved. Both ^dglap and Cdglap are known in first two orders of the perturbative QCD. 
Technically, it is simpler to calculate them at integer values of u> = n. In this case, the integrand of Eq. is called 
the n-th momentum of g^ s . Once the moments for different n are known, g NS at arbitrary values of u is obtained 
by interpolation. The expressions of the initial quark densities are obtained from phenomenological consideration, 
by fitting the experimental data at x ~ 1. Eq. @ shows that Cdglap govern the Q 2 - evolution whereas Cdglap 
evolve 8q{u>) in the x-space from x ~ 1 into the small x region. When in the ir-space the initial parton distributions 
8q(x) are regular in x, i.e. do not — > oo when x — > 0, the small- a; asymptotics of g± dglap is given by the well-known 
expression: 



9l DGLAP' 9l DGLAP ~ ex P 



Hl/x)ln(ln(Q 2 /^ 2 )/H^/A 2 QCD ) 



(3) 



On the contrary, when the total resummation of the double-logarithms (DL) and single- logarithms of x is doneQ, 
the Mellin representation for g^ s is 



g?»(x,Q*) = (e 2 /2) I —(l/xrCNs^Sq^expiHKsi^HQ 2 /^)) , 



(4) 



with new coefficient functions Cns, 



Cats(w) 



u>-Hfi>(u>) 



(5) 



and anomalous dimensions H^s, 



H NS = (1/2) u - y/uj 2 - B(u) 



(6) 



where 



B(w) = (4ttC f (1 +lj/2)A(u) + D{uj))/(2w 2 ) 



(7) 



D(ui) and A(uj) in Eq. Q are expressed in terms of p — ln(l/x), r] = ln(/i 2 /A.q CD ), b = (33 — 2nj)/127r and the color 
factors C F = 4/3, N = 3: 



D(u) 



2C F 



P + v i 
T — 

(p + rj) 2 + tt 2 r) 



(8) 



m = - 



dpe-^P 



r\ 2 + it 2 J (p + rj) 2 + tt 2 



(9) 
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H$ and Cns account for DL and SL contributions to all orders in a s . 
When x — > 0, 

g? s ~ (x 2 /Q 2 ) Ans/2 , gf ~ (z 2 /Q 2 ) As/2 (10) 

where the non-singlct and singlet intercepts are A^s = 0.42, As = 0.86. The x- behaviour of Eq. I|10|) is much steeper 
than the one of Eq. Obviously, the total resummation of logarithms of x leads to a faster growth of g\ when x 
is decreasing, compared to what is predicted by DGLAP, provided the input initial parton density Sq in Eq. (0) is a 
regular function of to at u> — > 0. 



III. ROLE OF THE INITIAL PARTON DENSITIES 



There are various forms in the literature for Sq(x), but all available fits include both a regular and a singular factor 
when x — > (see e.g. Refs. [l(J,lll| for detail). For example, one of the fits from Ref. is given by 



Sq(x) = Nr]x~ a (1 - xf '(1 + ^x") 



(11) 



with TV, rj being normalization factors, 
sum of pole contributions: 



a 



0.576, P = 2.67, 7 = 34.36 and 5 = 0.75. In the w -space Eq. (pH is a 



5q{uS) = Nrj (u — a) 1 + 2J m fc( w + Afe)~ 



(12) 



with Xk > 0, and the first term in Eq. l|12f) corresponds to the singular factor x~ a of Eq. (jlip. When Eq. I|ll|l is 
substituted in Eq. (|2J), the singular factor x~ a affects the small -x behavior of g\ and changes its asymptotics Eq. © 
for gi for the Regge asymptotics. Indeed, the small- x asymptotics is governed by the leading singularity lo = a, so 



gi DGLAP 



c(a)(i/xr((HQ 2 /A 2 ))/(Hv 2 /^)) 



7(a) 



(13) 



Obviously, the actual DGLAP asympotics of Eq. I|13f) is of the Regge type, and differs a lot from the conventional 
DGLAP asympotics of Eq. ©. Indeed it looks similar to our asymptotics given by Eq. (|10fl . namely by incorporating 
the singular factors into DGLAP initial parton densities it leads to the steep rise of g^ GLAP and therefore to a 
successful description of DGLAP at small x. In Ref. |9j it is shown that without the singular factor x~ a in the fit of 
Eq. 1)11(1 . DGLAP would not be able to work successfully at x < 0.05. In other words, the singular factors in DGLAP 
fits mimic the total resummation of logarithms of x of Eqs. Q4I10[I . To be more specific, although both (|13[1 and 1(1 OJI 
predict the Regge asymptotics for 171, there is a numerical difference in the intercepts: Eq. I)13|l predicts that the 
intercept of g± s should be a — 0.57, a value which is greater than our predicted non-singlet intercept Ajys = 0.42. 
Therefore the non-singlet g DGLAP grows, when x — > 0, faster than our predictions. Such a rise however is too steep 
and contradicts the results obtained in Refs. Q and confirmed in Refs. @. 



IV. CONCLUSIONS 



We have explicitly shown, by direct comparison of Eqs. (j3Jl and JT3J, that the singular factor x~ a in the Eq. I|ll|) 
for the initial quark density converts the exponential DGLAP-asympotics into the Regge one. On the other hand, 
comparison of Eqs. (|10|) and l|13(l also shows that this singular factor in the DGLAP fits mimics the total resummation 
of logarithms of x. This type of factors can be dropped when the total resummation of logarithms of x performed in 
Ref. [3 is taken into account. The remaining terms, which are regular in x in the DGLAP fits (the terms in squared 
brackets in Eq. I|ll[) ). can obviously be simplified or even dropped at small x and replaced by constants. A more 
detailed analysis as well as a suggestion to combine the leading logarithms resummation at small x with DGLAP 
can be found in Ref. |9j- The above results lead to an interesting conclusion: the expressions for the initial parton 
densities Sq used in DGLAP analyses have been commonly believed to be related to non-perturbative QCD effects; 
indeed they actually mimic the contributions of the perturbative QCD, so the whole impact of the non-perturbative 
QCD effects on g\ at small x is not large and can be approximated by a normalization constant. 
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